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ABSTRACT 

We present an analysis of the bulk composition of the MEarth transiting super Earth exoplanet 
GJ 1214b using planet interior structure models. We consider three possible origins for the gas layer on 
GJ 1214b: direct accretion of gas from the protoplanetary nebula, sublimation of ices, and outgassing 
from rocky material. Armed only with measurements of the planet mass (M p = 6.55 ± 0.98 M®), 
radius (R p = 2.678 ± 0.13 i?®), and stellar irradiation level, our main conclusion is that we cannot 
infer a unique composition. A diverse range of planet interiors fits the measured planet properties. 
Nonetheless, GJ 1214b's relatively low average density (p p = 1870±400 kgm -3 ) means that it almost 
certainly has a significant gas component. Our second major conclusion is that under most conditions 
we consider GJ 1214b would not have liquid water. Even if the outer envelope is predominantly 
sublimated water ice, the envelope will likely consist of a super-fluid layer sandwiched between vapor 
above and plasma (electrically conductive fluid) below at greater depths. In our models, a low intrinsic 
planet luminosity (< 2 TW) is needed for a water envelope on GJ 1214b to pass through the liquid 
phase. 

Subject headings: planets and satellites: general, planetary systems, stars: individual (GJ 1214) 



1. INTRODUCTION 

The era of super Earths is upon us with the first 
two transiting sub-Neptune mass exoplanets recently 
discovered. The first such tra nsiting planet, Co RoT- 
7b, is a M = 4.8 ± 0. 8 Mm (iLeger et all 120091 ) and 
R = 1.68 ± 0.09 i? ffi (|Queloz et all 120091) hot rela- 
tively dense planet with an average densit y similar to 
Earth 's. More recently, the MEarth project (jlrwin et al.l 
I2009D discovered transitin g low-mass planet GJ 1214b 
(|Charbonneau et al.l 12009). GJ 1214b has a mass of 
M p = 6.55 ± 0.98 M® and a radius of R p = 2.678 ± 
0.13 i?®. It is in a 1.5803952 ± 0.0000137 day pe- 
riod around an L« = 0.00328 ± 0.00045 L Q M dwarf 
of mass M* = 0.157 ± 0.019 M® and radius R* = 
0.2110 ±0.0097 R Q . 

GJ 1214b has a low enough density (p p = 1870 ± 
400 kgm -3 ) that it cannot be composed of rocky and 
iron material alone. The planet almost certainly contains 
a gas component. Even a planet of pure water ice is still 
too dense to match the observed mass and radius. At 
6.55 M® a pure zero-temperature water ice planet would 
have a radius of 2.29 i?® while an Earth-like composition 
would have a radius of about 1.64 i?®; these theoretical 
radii are 3 and 8<r lower than the value measured for 
GJ 1214b. While these simple arguments already reveal 
that GJ 1214b probably has a gaseous component, we are 
motivated to provide a more detailed analysis to quan- 
tify the range of possible planetary interior and gas layer 
compositions for GJ 1214b. 

We use planet interior structure models to constrain 
the bulk composition of GJ 1214b. In this work, we 
focus on three possible sources for the GJ 1214b gas 
layer: direct accretion of gas from the protoplanetary 
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nebula, sublimation of ices, and outgassing from rocky 
material. We examine end-member cases in which one of 
these three contributions dominates the gas layer. Based 
on GJ 1214b's mass and radius alone, we cannot infer 
a un i que interior composition (see, e.g . [ Valencia et al.l 
[20071 lAdanrs et al.l[200l iZeng fc Seagerll2008t ). Instead, 
there is a range of compositions that are consistent with 
the transit and radial velocity observations. Despite 
the inherent degeneracies plaguing the under-constrained 
problem of inferring GJ 1214b's composition from its 
mass and radius, we can nonetheless place interesting 
bounds on the gas envelope mass and draw insights into 
GJ 1214b's prospects for harboring liquid water. 

In Section[5]we explore the connection between the pri- 
mordial material comprising a planet and the sources of a 
planet's gas envelope. In Section[3]we describe our model 
of low-mass planet interiors. In Section|4]we present con- 
straints on the composition of GJ 1214b in each of three 
distinct scenarios for the origin of its gas layer. Discus- 
sion and conclusions follow in Sections [5] and [51 

2. CONNECTING GAS LAYER ORIGINS AND 
PLANET INTERIORS 

There is a wide range of possible chemical composi- 
tions for GJ 1214b's interior and gas layer. To motivate 
the discrete representative scenarios considered in this 
work, we look to the broad phases of materials that can 
contribute to a planet's bulk and to its gas layer. 

GJ 1214b may have formed from a variety of primor- 
dial material in its protoplanetary disk including gas 
(predominantly hydrogen and helium); ice-forming mate- 
rial (water, carbon monoxide, carbon dioxide, methane, 
and ammonia); and rocks or refractory material (iron, 
silicates, and sulfides). All three classes of primordial 
planet-building material (nebular gas, ices, and rocks) 
could have contributed to the gas layer observed on 
GJ 1214b today (Figure [J). Gas accreted directly from 
the nebula during planet formation, if retained, would 
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contribute hydrogen and helium. Sublimation of ices 
(for example, due to the release of gravitational energy 
during initial planet formation, the increased stellar ir- 
radiation following inward planetary migration, or late 
delivery of ices by comets) would produce H2O, CO, 
CO2, CH4, and NH3 vapor. Finally, rocky material 
can release volatiles to th e GJ 1214b gas layer via out- 
gassing during formation (jElkins-Tanton fc Seagerll2008t 
iSchaefer fe Feglevll2009bl) and tectonic activity after for- 
mation ( Kite et al.l 120091 ). Irrespective the origin of 
GJ 1214b's gas layer (be it from accreted nebular gas, 
sublimated ices, or outgassed rocky material), the gas 
envelope's mass and composition will have evolved over 
time under the influence of atmospheric escape. 

In this work, we focus on direct accretion of nebular 
gas, sublimation of ices, and outgassing from rocky ma- 
terial as possible sources for the gas layer on GJ 1214b. 
There are, however, other atmosphere formation pro- 
cesses worth mentioning. Vaporization of rocky material 
can contribute to atmospheres surro unding highly irra- 
diated super Earths like CoRoT-7b (jSchaefer fc F cglcv 
l2009aj) . but GJ 1214b is not hot enough for this process 
to occur. Even at temperatures too low to sublimate 
ices or vaporize oceans, volatiles stored as clathrate hy- 
drates in icy material can be outgassed into a planet's 
atmosphere. Sputtering by the stellar wind and microm- 
eteorites, photolysis, radiolysis, and chemical reactions 
between stellar wind ions and planet surfaces all con- 
tribute to tenuous atmospheres surrounding solar system 
bodies. The contributions of these gas sources are neg- 
ligible, however, compared to the gas volume needed to 
account for GJ 1214b's transit depth. 

The actual bulk make-up of GJ 1214b is determined 
by its unknown formation, migration, and evolution his- 
tory. For instance, if GJ 1214b initially formed beyond 
the snow line it would contain more icy material than 
if it formed closer to its star. The mass of nebular 
gas initially captured would depend upon the nascent 
GJ 1214b's accretion luminosity as well as the local con- 
ditions (density, temperature, opacity, and mean mol ecu- 
lar weight) in the protoplanetary disk {Rafikov 2003). In 
the solar system, there is a definite relationship between 
the relative abundances of rock-ice-gas and planet mass: 
small planets (< 1 M s ) are rocky, intermediate planets 
(~ 15 — 17 M©) are icy, and larger planets are predom- 
inantly composed of H and He. Ro ugh constraint s on 
solar system planet compositions from lGuillotl (|2005| ) are 
plotted in Figure [2j We do not attempt to tighten the 
constraints on GJ 1214b's interior by directly incorporat- 
ing planet formation theories into our analysis. Instead 
we allow for the full range of primordial gas-ice-rock ra- 
tios and explore the constraints imposed on these ratios 
by the measured mass and radius. 

We consider a series of scenarios that encompasses all 
nebular gas-ice-rock mass fraction combinations for the 
primordial material making up GJ 1214b (Figure [5]). In 
each case, we assume that a single contributor (nebular 
gas, ice, and rock) dominated as the source for the gas 
envelope observed today on GJ 1214b. We are thus con- 
sidering end-member scenarios within the continuum of 
possible gas envelope compositions. In case I, GJ 1214b 
managed to accrete and retain H and He gases from the 
nebula, and includes primordial gas, ice, and rock in its 
bulk make-up. In this scenario, we neglect any contribu- 
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Fig. 1. — Diagram detailing possible sources considered in this 
work for GJ 1214b's gas layer. The segments of the circle rep- 
resent the three categories of primordial material that could have 
contributed to forming GJ 1214b: refractory materials, ice-forming 
material, and nebular gas. All three categories of primordial ma- 
terial can contribute to an eventual planetary gas layer. The gas 
formation processes we consider for GJ 1214b are indicated by the 
black arrows. Each arrow points to a box describing possible ini- 
tial chemical compositions for the gas contributed by e ach source. 
The ice compositions were taken from Marbocuf et al. (2008), and 
the outgassed at mosphere compositions refl ect the chemical equi- 
librium results of Schacfcr & Fcglcy (2009b). The gas layer com- 
position will evolve over time under the influence of atmospheric 
escape. While primordial gas, ice, and rocks are all given equal 
fractions of the circle in this diagram, the primordial gas-ice-rock 
ratios of planets can vary over a wide range and affect the relative 
importance of the three gas layer sources shown in the diagram. 
In this work, we consider three end-member cases (labeled as cases 
I, II, and III) in which a single gas layer source dominated on 
GJ 1214b; the case associated with each gas formation process is 
indicated in the diagram. 
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Fig. 2. — Schematic diagram illustrating the range of possible 
planet primordial bulk compositions. In this figure "gas" refers to 
primordial H and He accreted from the nebula, "ice" refers to ice- 
forming materials, and "rock" refers to refractory materials (e.g., 
iron and silicates). Constraints on the current compositions of the 
solar system planets are plotted in purple (planets are denoted 
by their first initial). For GJ 1214b, we consider the full possible 
range of primordial gas, ice and rock relative abundances. In case 
I (green), GJ 1214b accreted and retained primordial gas, ices and 
refractory material. In case II (blue), GJ 1214b did not retain any 
primordial gas. incorporating only icy and rocky materials. Finally 
in case III (red), GJ 1214b formed from purely rocky material. 
Planets in cases II and III, which do not contain any primordial 
gas, may still harbor a gas layer produced b y sublimated ices or by 
outgassing . This diagram was inspired by Chambers (201(3) and 
Ste venson (2004). 



3 



tions to the gas envelope from the ices or rock, and take 
the current gas envelope on GJ 1214b to be composed of 
H and He. In case II, GJ 1214b has an interior formed 
from icy material and rock and did not retain any nebular 
gas (either having never accreted any in the first place or 
having lost any nebular gas that it once had). Here, we 
assume that vapors from the icy material dominate the 
gas layer (neglecting any contributions from the rocky 
material). Finally, in case III, GJ 1214b formed from 
purely rocky material and did not acquire any ices or gas 
from the protoplanetary disk. In the absence of accreted 
ices and gas, the planetary gas envelope must originate 
by outgassing during formation or tectonic activity. 

3. STRUCTURE MODEL 

Our model for the interior st ructure of low-mass exo - 
plancts is described in detail in lRogers fc Seagerl (|2010f ). 
Here, we give a brief summary. 

We assume spherically symmetric and differentiated 
planets consisting of a core, mantle, ice layer, and gas 
envelope. The equation for the mass of a spherical shell 
and the equation describing hydrostatic equilibrium form 
a coupled set of differential equations for the radius r (to) 
and pressure P(m), viewed as functions of the interior 
mass m, 



dr 
dm 
dP 
dm 
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where p is the density and G is the gravitational constant. 
The equations of state (EOS) 

Pi = fi(P,T) (3) 

relates the density p (m) to the pressure P (m) and tem- 
perature T (to) within each distinct chemical layer i. 
We integrate Equations [1] and [2] imposing that both 
P and r axe continuous across layer boundaries. The 
outer boundary conditions on the pressure and opti- 
cal depth are calculated fo l lowing the procedure de- 
scribed in IRogers fc Seagerl (]2010l ) so as to take int o 
account the 'transit radius effect' (jBaraffe et al.l 12003). 
We then solve iteratively for the core-mantle mass ra- 
tio that yields a consistent solution (for a given distri- 
bution of mass in the outer layers, set of atmospheric 
parameters, total planet mass and radius). This is 
the common general approach fo r modeling planet in- 
teriors (e.g. iValencia et all 120061: iFortney et all 120071: 



Adams et all 120081 [Zeng fc Seagerl 120081: Baraffe et al.l 



20081: iFigueira et all 120091: iGrasset etafl 120091 ). The 
model we use in t his work is improved over 
IRogers fc Seagerl (|2010l ) by including a temperature- 
dependent water EOS and is different from most pre- 
vious planet interior models by providing quantitative 
constraints on the range of plausible gas envelope masses 
for a given planet mass and radius. 

The thermal profile of our model planets is divided into 
three regimes: an outer radiative regime in the gas/fluid 
envelope, an inner convective regime in the gas/fluid en- 
velope, and a solid interior in which thermal effects are 
neglected. We assume that in the outermost region the 
planets' gas/fluid envelopes are in radiative equilibrium, 
and use a temperature profile derived from an analytic 



"two-stream" solution to the gray equations of radia- 
tive transfer for a plane-parallel irradia ted atmosphere 
(Equation (45) in iHansenl |200g) . The |Hansenl (|2008h 
temperature profile describes the temperature T (r) as 
a function of the optical depth, r, and depends upon 
the degree of stellar insolation, internal luminosity of the 
planet, and the ratio of the thermal to visible opacities 
(parameterized by To, T e ff, and 7, respectively). The on- 
set of convective instabilities (0 < (dp/ds) p ds/dm) de- 
limits the transition to the convective layer of the fluid 
envelope. In the convective regime, we adopt an adia- 
batic temperature profile. Finally, within the solid ice, 
mantle, and core of the planets we neglect the tempera- 
ture dependence of the EOSs, employing an isothermal 
temperature profile. At the high pressures found in the 
solid interior layers, therma l corrections have o nly a small 
effect on the mass-density (Scager et al. 2003). 

We choose both the fiducial values and uncertainty 
ranges for the atmospheric parameters in our model 
(Tp, T e g, and 7) follow ing the prescription described in 
IRogers fc Seagerl (|2010l ). As our fiducial value of To, we 
take 558 K, the equilibrium temperature of GJ 1214b as- 
suming full redistribution and neglecting reflection. We 
also consider a range of To values from 789 to 501 K, re- 
flecting uncertainties in the planet's albedo and in the de- 
gree to which energy is redistributed within the planet's 
gas envelope. Planetary bond albedo values up to 0.35 
are considered. The parameter T c g describes the intrinsic 
luminosity of the planet. GJ 1214's old-disk kinematics 
and lack of chromosphcric activity suggest that it has 
a ste llar age between 3 and 10 Gyr (|Charbonneau et all 
2009) . We employ a simple approximate scaling relation, 
derived from a powe r-law fit to cooling calculations from 
IBaraffe et al.l (|2008l ), to relate the planets intrinsic lumi- 
nosity to its mass, radius, and age. For a planetary age of 
3-10 Gyr, we estimate a plausible range for the intrinsic 
luminosity of the planet T e g = 44 — 69 K, while adopting 
a fiducial value of T c g = 61 K corresponding to 4.5 Gyr. 
For gas envelopes composed of some mixture of H and 
He, we adopt a fiducial value of 7 = 1 and consider a 
range from 7 = 0.1 to 10. Because the opacity of water 
is far higher in thermal wavelengths than in the visible, 
we expect that 7 < 1 in a water envelope. For our water 
vapor atmospheres we thus adopt a fiducial value 7 = 0.1 
and consider an uncertainty range from 7 = 0.01 to 1. 

In the solid (uniform temperature) l ayers of the 
planet we employ EOS data sets from ISeager et al.l 
(2007) which were derived by combining experimental 
data at P < 200 GPa with the theoretical Thomas- 
Fermi-Dirac equation of state a t high pressures ; P > 
10 4 GPa. We consider Fe (e) (I Anderson et al.l |200C 
Mg 1 _ x Fe x Si0 3 perovskite (Elkins- Tantonl I2008D . and 
H 2 ice (jHemlev et all 119871) . To describe hydro- 
gen and heliu m envelopes we use EOSs for H/He 
mixtures f rom ISaumon et all (119951 ) a nd opacity ta- 
bles f rom iFreedman et al.l ~ 2008[ ) and iFerguson et al.l 
(2005). We have compiled a temperature-dependent 
EOS for water up to 32 GPa spanning liquid, vapor, 
super-fluid, and plasma phases. Our water EOS com- 
bines data from the "The IAPWS Formulation 1995 
for the Thermodynamic Properties of Ordinary Water 
Substance for Genera l and Scientific Use" (IAPWS-95; 
iWagner fc Prufj|2002l) retrieved from the NIST Chem- 



4 



Rogers et al. 



istry WebBook (|Lemmon et al.l l2OO90. extrapolations of 
the IAP WS-95 formulat ion calculated using FLUIDCAL 
software (Wagner 2009 1), and th e Thomas-Fermi-Dirac 
EOS (Salpetcr & Zapolsky 1967). For the opacity in the 
water vapor layer we use Planck means calcula ted with 
molecular line data from iFreedman et al.1 ([2008). 

4. RESULTS 

4.1. Background 

We consider three cases for the interior makeup of 
GJ 1214b: I) a planet that formed from nebular gas, ice, 
and rock and still harbors a primordial H/He envelope; 

II) an ice-rock planet that failed to accrete H/He gas from 
the protoplanetary disk but now has a vapor envelope; 

III) a rocky planet with an outgassed atmosphere (Fig- 
ures [1] and [2]) . These scenarios determine what distinct 
chemical layers we consider in our differentiated planet 
model described in Section[3] We assume Y = 0.28 H/He 
for the nebular gas, and pure H2O ice for the icy mate- 
rial. We model the rocky material by a combination of 
metallic iron and Mgo.gFeo.iSiOa silicates without im- 
posing any a priori constraints on the iron-to-silicates 
ratio. We assume that, during GJ 1214b's formation, 
the primordial rocky material differentiated to form an 
iron core and silicate mantle in the planet. In this way, 
the rocky material contributes two layers in our planet 
interior structure model. The effect of choosing other 
chemical compositions to represent the primordial gas, 
ice, and rock is discussed in Section 15.51 

We say that an interior composition is consistent with 
the measured planetary mass and radius within their 
na observational uncertainties if there is some choice 
of planet mass within (M p — no~M p , M p + ugm v ) , planet 
radius within \R V — no~n v , R p + no~R p J , and atmospheric 
parameters (T c g, To, and 7) within the ranges given in 
Section [3] that yields a consistent solution for the interior 
composition. A mor e sophisticated error an alysis (such 
as that described in IRogers fc Seagerl I2010D is not yet 
warranted given the current error bars on GJ 1214b's 
mass and radius. 

We employ ternary diagrams to graphically present our 
composition constraints for each scenario (Figures |3] and 
3]). Ternary diagrams are useful tools to graphically rep- 
resent three component data (x, y, z) for which the com- 
ponents are constrained to be positive (x, y,z > 0) and 
to have a constant sum (x + y + z = 1). Since such data 
have only 2 degrees of freedom, it could easily be dis- 
played with a x-y Cartesian plot, wherein the axes and 
the z — line would form a right triangle. To show all 
three components (x, y, z) on an equal footing, the x-y 
Cartesian plot can be squished (via a linear transforma- 
tion) so that the x and y axes meet at a 60° angle and 
form an equilateral triangle with the z = line. The 
resulting equilateral triangle diagram is a ternary dia- 
gram. The three vertices of the diagram represent points 
where x = 1, y = 1, and z = 1, while x — 0, y = 0, and 
z = along the respective opposing edge. At each inte- 
rior point, the value of x is given by the perpendicular 
distance from the x = edge, with the values of y and 
z defined analogously. More detailed desc riptions of how 
to rea d ter nary diagrams can be found in Valenci a et al.l 
(f2007h and lZeng fc Seagerl (|200l) . 



4.2. Case I: Gas-Ice-Rock Planet with Primordial Gas 

Envelope 

We first consider the case in which GJ 1214b managed 
to acquire and retain H/He gas from the protoplanetary 
disk. In this scenario GJ 1214b incorporated primordial 
iron, silicates, ice, and gas into its bulk make-up. We 
thus allow for four chemically distinct layers in the planet 
interior: an iron core, silicate mantle, water-ice layer, and 
H/He envelope (with Y — 0.28). 

Despite the range of allowed compositions, the mass 
of GJ 1214b's H/He gas envelope is tightly constrained 
(Figure [3]) . A gas mass fraction of 0% is not allowed 
within the la observational error bars on M p and R p . 
The maximum and minimum gas envelope masses in this 
scenario occur for end-member mass distributions where 
the gas envelope surrounds a pure iron interior and a 
pure water ice interior, respectively. For M p = 6.55 M e , 
R p = 2.678 R®, the mass in the gas layer could be at 
most 3.6%-5.2% of the planet mass depending on the 
atmospheric thermal profile (hotter atmospheres require 
less mass of H and He to occupy a similar volume) . When 
M p and R p are varied within their la observational un- 
certainties, the range of maximum H/He mass fractions 
widens to 3.2%-6.8%. The minimum gas mass fraction is 
more sensitive than the maximum to the atmospheric en- 
ergy budget and composition (e.g., metallicity and opac- 
ities) and could be in the range 9 x 10 -5 and 2 x 10~ 3 
at the fiducial planetary mass and radius. 

In this scenario, GJ 1214b has a less massive H/He 
envelope than our solar system Neptune, whose compo- 
sition is roughly 5% -15% H and He, 60%-70% ices, and 
25% rocks by mass (jPodolak et alj|1991t iHubbard et all 
1995). GJ 1214b could, nonetheless, support a gas en- 
velope that is large as compared to the terrestrial so- 
lar system planets. For a Ganymede- like interior with 
iron:silicates:water ice in the ratio 3:22:75 by mass, 
GJ 1214b requires an H/He envelope accounting for be- 
tween 0.01% and 0.6% of the planetary mass depend- 
ing on the atmospheric temperature. This is up to 60 
times larger than the atmosphere mass fraction on Venus 
(~ 10- 4 ). 

4.3. Case II: Ice-Rock planet with Sublimated Vapor 
Envelope 

A planet interior dominated by ice and a concomitant 
gas envelope dominated by vapors from ice-forming ma- 
terials is an intriguing possibility for GJ 1214b. This 
scenario is substantially different from any of the solar 
system planets, but could be thought of as a class of 
bigger , hot ter versions of Jup iter's icy moons. iKuchnerl 
(2003) and lLeger et al.l (|2004D first proposed that water- 
rich planets might be prevalent on orbits accessible to 
transit and radial velocity detections, although no such 
planets have been conclusively discovered so far. The for- 
mation pathway for these planets involves inward migra- 
tion of proto-planets that formed from volatile ice-rich 
material beyond the snow line but that never attained 
masses sufficient to accrete large amounts of H/He neb- 
ular gas. If GJ 1214b falls into this category, it could be 
the first discovered member of a whole new population 
of exoplanets. 

For this scenario, in which GJ 1214b did not accrete 
or retain any H or He from the protoplanetary disk, we 
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Fig. 3. — Ternary diagram for case I in which GJ 1214b re- 
tained a primordial H/He envelope, having formed from primor- 
dial gas, ice, and rock. The relative contributions of the iron core, 
Mgo.9Feo.iSi03 mantle, H2O ices, and H/He envelope to the mass 
of the planet are plotted. The core and mantle are combined to- 
gether on a single axis, with the vertical distance from the upper 
vertex determined by the fraction of the planet's mass in the two 
innermost planet layers. The black shaded region denotes the inte- 
rior compositions that are consistent with the nominal planet mass 
and radius (Mp = 6.55 M®, R p = 2.678 -R©) for our fiducial choice 
of atmospheric parameters (7 = 1, To = 558 K, T e jj = 61 K). The 
H/He mass fraction has a spread in this case due to the range of 
possible core-to-mantle mass ratios. The span of plausible interior 
compositions widens to the green shaded area when the range of 
atmospheric parameter values delimited in Section [3] is considered. 
The red, yellow, and blue shaded regions denote compositions that 
are consistent with M p and R p to within 1, 2, and 3<r of their 
observational uncertainties, respectively, when uncertainties in the 
atmospheric parameters are also included. 

adopt an interior planet structure consisting of an iron 
core, silicate mantle, and water envelope. The pressure- 
temperature (PT) profile determines the phase of water 
in the envelope; vapor, liquid, super-fluid, high pressure 
ices, and plasma phases are all included in our H2O EOS. 
We model the thermal profile of the water envelope fol- 
lowing the same prescription as for the H/He layers (de- 
scribed in Section [3J . To allow for the presence of a 
greenhouse effect, we take 7 = 0.01 — 1 for H2O (com- 
pared to 7 = O.f — 1.0 for H and He). 

We show in Figure 0] the possible distributions of 
mass between the core, mantle, and water envelope that 
are consistent with the measured mass and radius of 
GJ 1214b. A sublimated vapor dominated envelope on 
GJ 1214b is possible if water accounts for a large frac- 
tion of planet mass. At the fiducial measured planet 
mass and radius (M p = 6.55 M®, R p = 2.678 R®) at 
least 88% H 2 by mass is required. To account for the 
observed planet mass and radius within their 1, 2, and 
3cr observational uncertainties, at least 47 %, 24%, and 
6% wa ter by mass are require d, respectively. ILeger et al.l 
(2004) and Selsis et al. (2007) have put forward that ice- 
rock planets formed beyond the snow line may have a 
comet-like bulk composition with 50% H2O and 50% sil- 
icates and iron by mass. This composition is consistent 
with the measured GJ 1214b radius within la. Interior 
structure considerations do not preclude the possibility 



that GJ 1214b is water rich. 

GJ 1214b does not contain liquid water in any of our 
model interiors displayed in Figure |4j The PT profiles 
that result from the range of equilibrium temperature 
and internal heat flux (To and T e g) values we consid- 
ered are too hot to allow liquid water, even at high pres- 
sures. Our putative GJ 1214b water envelopes begin 
in the vapor phase at low pressures, then continuously 
transition to a super-fluid at P = 22.1 MPa (the crit- 
ical pressure of water), before eventually becoming an 
electronically conductive dense fluid plasma at greater 
depths (T > 4000 K). 

To obtain liquid water in our model interior, we must 
decrease GJ 1214b's assumed intrinsic energy flux below 
5 x 10 -4 W m~ 2 or T er j < 10 K; fo r comparison, Earth's 
internal heat flux is 0.087 W m" 2 (|Turcotte fc Schubert! 
2002;). The intrinsic luminosity of GJ 1214b is very un- 
certain, and a detailed evolution calculation to better 
constrain its magnitude is out of the scope of this work. 
Nonetheless, we predict that, in this scenario, GJ 1214b 
would need a cold interior in order to harbor liquid water. 

Our conclusions regarding the possibility of liquid 
water in an extended vapor atmosphere on GJ 1214b 
are contingent upon our parameterized PT profile ade- 
quately describing the water envelope. It is important to 
note that our model thermal profile assumes GJ 1214b is 
in radiative equilibrium with the incident stellar irradia- 
tion it receives at its current orbital location. If GJ 1214b 
has recently migrated and is undergoing active vaporiza - 
tio n (as in the scenarios considered by iKuchnerl (2003) 
and lValencia et all (|2009f )) it could have temperature and 
H20-phase profiles within its envelope that are very dif- 
ferent from those available within the framework of our 
models (e.g., atmospheres with a temperature inversion 
or those out of radiative-convective equilibrium). 

4.4. Case III: Rocky Planet with Outgassed 
Atmosphere 

We turn to the possibility that GJ 1214b formed from 
purely rocky material without retaining any H/He gas 
or icy material from the protoplanetary disk. We reiter- 
ate that GJ 1214b must still have a substantial gas layer 
in this case, because a rocky planet is too dense to ac- 
count for the measured mass and radius (e.g., a gasless 
Earth-like composition yields a planet radius that is 8a 
lower than that measured for GJ 1214b). In this rocky 
planet scenario, an outgassed atmosphere contributes to 
the GJ 1214b transit radius. We focus here on an atmo- 
sphere produced by outgassing during planet formation; 
outgassing from post-formation geological activity is an- 
other possibility, which is discussed in Section 15.31 

It is difficult to predict a priori the composition of 
the gas layer that would be produced by outgassing 
on GJ 1214b. The initial composition of the out- 
gassed atmosphere is strongly dependent on the com- 
position of planetcsimals comprising GJ 1214b. Even 
among solar system chondrites, the outgassed atmo- 
sphere compositions would range from ^-dominated, to 
H2Q-dominated, to CO-do minated, to C02-dominated 
(|Schaefer fc Feglevl I2009bf ). In addition, the condi- 
tions during magma solidi fication and outgassing affect 
the gas layer outcome (Elkins-Tanton fc Seagerl 120081 : 
iSchaefer fc Feglevl 12009b ). Further complicating the 
rocky planet gas layer, is the subsequent atmospheric 
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Fig. 4. — Ternary diagram for case II in which GJ 1214b formed 
from refractory material and ices and has an envelope dominated 
by vapor from sublimated ice. The fractions of the planet's mass 
in the Fe core, the Mg 9 Feo.iSi03 silicate mantle, and the water 
vapor envelope are plotted on the three axes. The solid black curve 
represents the locus of interior compositions that are consistent 
with the nominal planetary mass and radius (M p = 6.55 Mj, 
R p = 2.678 i?e) for our fiducial choice of atmospheric parameters 
( 7 = 1, T = 558 K, T cff = 61 K). The colors in this figure have 
the same designations as in Figure [3] 

escape (|Charbonneau et all I2009H and photochemistry. 
Despite this uncertainty in composition, we can never- 
theless make some concrete statements about a putative 
outgassed envelope on GJ 1214b. 

If GJ 1214b's gas layer was produced by outgassing, a 
substantial fraction of it must be in a component that 
is less dense than water vapor. The light component is 
needed for GJ 1214b's envelope to have both a mass low 
enough to be produced by outgassing and a volume large 
enough to account for the transit radius. From our study 
of case II, we found that if water is the least dense compo- 
nent of GJ 1214b, at least 47% H2O by mass is required 
within la (Figure 2]). This is more water than a terres- 
trial planet formed from chondritic planetesimals can de - 
gas (up to 23% by mass) (|Elkins-Tanton fe Sea"geill2008f) . 
Volatile molecules heavier than water (such as CO, CO2, 
and N2) have even smaller upper bounds (in % planet 
mass) on the amounts they can be outgassed. The re- 
quirement for substantial quantities of a light species is 
interesting because it limits the outgassed atmosphere 
compositions that can be relevant to GJ 1214b from 
among the wide range of a priori possibilities. In partic- 
ular, an Earth-like ^-dominated outgassed atmospheric 
composition and a Venus-like CO2 outgassed atmosphere 
are both ruled out for GJ 1214b. 

Molecular hydrogen, H2, is the most likely can- 
didate for an atmospheric species that is both 
light enough and plausibly outgassed in sufficient 
quantities to account for GJ 1214b's transit ra- 
dius. In fact, H2 is predicted to dominate the 
atmospheres outgassed by ordinary H, L, LL, and 
high ir on enstatite EH chondrite-composition planetes- 
imals (jElkins-Tanton fe Seagerl [20081: iSchaefer fe Feglevi 
l2009bl) . Although He has a low molecular weight, it 



does not bind to minerals the way H does, and con- 
sequently it cannot be accreted with the rocky pri- 
mordial material and later rele ased through outgassing 
(jElkins-Tanton fe Seagerl 120081 ). The other possible 
species, CH4 and NH3, have molecular weights only 
slightly lower than water and are not typically outgassed 
in large quantities. 

We emphasize that only a relatively small "out- 
gassable" amount of H2 is required to account for 
GJ 1214b's transit radius. We show this quantitatively 
by considering a pure hydrogen gas layer surrounding a 
rocky interior in Figure [5] Just 5 x 10~ 4 of the planet 
mass in H2 surrounding a rocky core is sufficient to 
account for the transit radius to within la. This is 
2 orders of magnitude below the maximal amount of 
hydrogen that might be outgassed by a rocky planet 
(6% of the planet mass), which corresponds to the ex- 
treme where the planet formed from planetesimals simi- 
lar to EH chondrites and fully o xidized during formation 
( Elki ns-Tanton fe Seagerl [2001 . For an iron core mass 
fraction similar to Earth's (30%), GJ 1214b needs 0.3%- 
1.2% of its mass in a pure H2 gas layer at the nominal 
mass and radius. This H2 gas layer is small compared 
to the H2O envelopes in the water planet scenario, but 
is still large compared to Earth's atmosphere (which ac- 
counts for roughly 0.0001% of Earth's mass). 

A pure H2 envelope as we have assumed above is some- 
what artificial; realistically H2 will not be outgassed on 
its own but in combination with heavier molecules. At 
a given Fe core mass fraction, including H2O, CO, CO2, 
and other additional species in the outgassed envelope 
would in general tend to decrease the mass fraction of 
H2 required to reproduce the planet radius (Figure [5J , 
while increasing the total gas mass fraction of all volatile 
species combined. This is because i) the heavier out- 
gassed volatiles are still less dense than silicates and ii) 
the additional atmospheric species would increase the 
opacity of the atmosphere. This strengthens the result 
that an outgassed hydrogen-rich envelope surrounding a 
rocky solid interior is not ruled out by the measured mass 
and radius of GJ 1214b. We consider how this interpreta- 
tion is influenced when atmospheric escape is taken into 
account in Section I5T2"1 

5. DISCUSSION 
5.1. Atmosphere Observations 

What further observations can discriminate among the 
gas layer origin possibilities? A promising possibility is 
spectral observations of the planetary atmosphere. Here, 
we summarize some ideas for how atmospheric spectra 
might help to distinguish between a gas envelope domi- 
nated by nebular gas, sublimated ices, or outgassing. 

Atmospheric He is the discriminator between 
captured H/He envelo p es o r outgassed envelopes 
(|Elkins-Tanton fe Seagerl 120081) . He is difficult to 
observe spectroscopically at planetary temperatures, 
however. Recent outgassing by geological activity could 
be revealed through the detection of spectral features 
from species with sho rt photochemical lifetimes such as 
S0 2 (IKite et al.H2009l) . A water vapor dominated atmo- 
sphere would have saturated water vapor features (e.g., 
absorption b ands a t A ~ 5 — 8/xm and A ~ 16 — 18/im 
iLeger et"al"1 ((2004)), but these may not be unique 
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Fig. 5. — GJ 1214b as a rocky planet with an outgassed atmo- 
sphere. A three-layered planet structure composed of an iron core, 
silicate mantle, and pure hydrogen envelope is assumed. The frac- 
tion of the planet's mass in the hydrogen gas envelope is plotted 
as a function of the fraction of the planet's mass in the iron core. 
Any mass not in the hydrogen envelope or iron core is contained 
in the silicate mantle; the dashed black lines represent contours of 
constant mantle mass fraction. The colored shaded regions have 
the same designations as in Figure [3] 

identifiers. Discriminating between a hydrogen-rich 
envelope and a water vapor at mosphere on the basis o f 
the atmospheric scale height ( Miller- Ri cci et all I2009T) 
might be more promising than discriminating on the 
basis of the presence or absence of spectral lines. The 
scale height is accessible via the depth of strong (or sat- 
urated) transit transmission spectral features. Although 
it will be tricky to constrain interior compositions from 
atmospheric spectra, the prospect of greater insights 
into the composition and formation of GJ 1214b and its 
envelope is compelling. 

5.2. Atmospheric Escape 

We cannot definitely rule out any of our three compo- 
sition scenarios on the basis of atmospheric escape cal- 
culations. The water vapor envelope in case II is most 
resilient against atmospheric escape. While mass loss 
is a more important consideration for the hydrogen-rich 
primordial and outgassed envelopes in cases I and III, 
uncertainties in the mass-loss rate and planet age leave 
room for the possibility that GJ 1214b retained sufficient 
H/He or H2 to account for its transit depth. Further ob- 
servations are warranted to tighten the constraints on 
the origin and character of GJ 1214b's gas layer. 

Atmospheric escape may have an important influence 
on mass and composition evolution of the GJ 1214b 
gas envelope. We calculate 2 x 10 6 kgs -1 for an 
energy-limited upper bound on the current mass-loss 
ra te of GJ 1214b follow i ng th e approximate approach 
of lLecavelier Pes Etangsl (|2007t ). This is on the same or- 
der as the 9 x 10 5 kgs -1 escape rate iCharbonneau et al.l 
(2009) predicted for a hydrogen-rich atmosphere escap- 
ing hydrodynamically from GJ 1214b. At 2 x 10 6 kgs -1 , 
0.03 M® would be lost over 3 Gyr and 0.1 Mq would be 
lost over 10 Gyr. The actual cumulative mass lost over 



GJ 1214b's lifetime may be much higher than these val- 
ues because the host star was probably brighter in UV 
at earlier times. 

Out of the three cases we considered for GJ 1214b's 
composition, case II, the ice-rock scenario having a water 
vapor envelope, is most robust against atmospheric es- 
cape. First, water has a higher molecular weight than H 
or He which makes it easier to retain (although photodis- 
sociation of H2O may be important). Second, with H2O 
comprising upward of 50% of GJ 1214b's current mass 
in this scenario, a cumulative water loss on the order of 
0.1 will not have significantly changed the overall 
character of the planet. In contrast, losses of this mag- 
nitude are far more significant for the hydrogen-rich en- 
velopes considered in cases I (primordial gas dominated) 
and III (outgassing dominated) , since these envelopes can 
account for at most a few percent of GJ 1214b's current 
mass. 

In case I, atmospheric escape increases the amount of 
gas GJ 1214b would have needed to accrete from the pro- 
toplanetary nebula. Protoplanetary cores between 1 and 
10 M ffi can acquire substantial primordial atmospheres 
even if the cores are too small for the nucleate d insta- 
bility and runaway gas accretion t o commence (|Mizunol 
1980). For instance, iRafikovi ((2006h found that a 6.5 M e 
core in a minimum mass solar nebula could develop an 
atmosphere of up to several tenths of an Earth mass of 
nebular gas (depending on the core's semimajor axis and 
accretion luminosity). Whether GJ 1214b could accrete 
and retain enough nebular gas to account for its transit 
radius is ambiguous. 

Atmospheric mass loss creates the most tension with 
case III, because there is an upper bound on the cu- 
mulative amount of gas GJ 1214b can outgas over its 
lifetime. Based on solar system meteorites, a rocky 
planet can outgas at most 6% of its mass as hydrogen 
( Elki nsTTanton fe Seagerl2008l) and less will be outgassed 
if the iron content of the rocky material is not fully oxi- 
dized. It is likely that GJ 1214b has lost most of any Pi- 
atmosphere outgassed during formation. However, only 
~ 0.1% of the planet mass in H2 surrounding a rocky 
core could be needed to account for the transit radius to 
within la. Due to the large uncertainties on the planet 
age and time averaged mass-loss rate, the outgassing can- 
not be fully ruled out as the source of GJ 1214b's gas 
layer. 

5.3. Ongoing Outgassing 

If GJ 1214b is a rocky planet with an outgassed at- 
mosphere as in case III, ongoing outgassing by geolog- 
ical activity could be contributing to its gas envelope. 
However, the rate of ongoing outgassing is e xpected to 
be sm aller than the atmospheric escape rate. Kite et al~l 
(|2009f ) predict that a - 6.5 M e super Earth experienc- 
ing plate tectonics would have rates of volcanism per 
unit planet mass 6 times Earth's current rate at planet 
age of 3 Gyr and 0.2-0.3 times Earth's current rate at 
a planet age of 10 Gyr. Lower rates of volcanism are 
predicted for planets 3-10 Gyr old that are not tectoni- 
cally active. If we consider the most optimistic case for 
the rate of volcanism of GJ 1214b and assume a magma 
volatile content similar to volatile-rich terrestrial mag- 
mas at mid-ocean ridges (1.5 % H2O and 400 ppm CO2 
by mass lOppenheimerl (|2003D ). we obtain an upper limit 
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on the volcanic outgassing rate of 1 x 10 6 kgs -1 H2O and 
4 x 10 4 kgs -1 CO2. Although outgassing by geological 
activity may help to replenish the GJ 1214b envelope, 
it probably cannot completely offset the effect of atmo- 
spheric escape if the planet is older than 3 Gyr. 

5.4. Necessity of a Gas Layer 

GJ 1214b's low average density implies that it has a 
low density gas envelope. This statement is valid as 
long as the true planet mass and radius lie within the 
2a measurement uncertainties. If GJ 1214b's mass and 
radius both differ from their measured values by more 
than ~ 2<T, a solid, ice-dominated interior composition 
with no gas envelope is barely allowed. A pure ice 
planet seems physically implausible, however, because 
silicates (i.e., higher density material) are expected to 
be accreted along with i ces during planet formation. 
iCharbonneau et "all (|2009f ) point out, however, that the 
stellar radius they derived for GJ 1214 from observations 
is 1 5% larger than that predicted by the theoretical mod- 
els IBaraffe et al.l ([T998). If systematics have led to an 
overestimation of the planet radius, the evidence for a 
gas layer on GJ 1214b would be reduced. 

5.5. Model Uncertainties 

Despite our quantitative constraints on the range of 
interior compositions, we are faced with some uncertain- 
ties. One uncertainty is the atmospheric temperature, 
controlled by the unknown interior energy and unknown 
albedo. A hotter atmosphere fills a larger volume than 
a cooler atmosphere, requiring less atmospheric mass to 
fit the planet radius. We have chosen a range of rea- 
sonable values for the parameters governing the atmo- 
spheric PT profile (Section [3]) and found that the model 
uncertainty is roughly comparable to the observational 
uncertainties. Although we adopted nearly identical at- 
mospheric parameter ranges for all three interior struc- 
ture cases we considered, in reality GJ 1214b's albedo, 
atmospheric absorption, and interior luminosity are cou- 
pled to its interior composition. 

The precise chemical make-up of GJ 1214b's interior 
layers is another source of uncertainty. We have adopted 
artificially clean single-chemical materials to represent 
GJ 1214b's core (pure Fe), mantle (Mgo.gFe .iSi03 per- 
ovskite), and ice layer (pure H2O). In reality, we expect 
a mixture of chemical compounds to contribute to each 
layer. The presence of a light element (such as sulfur) 
in the iron core would decrease the mass required in the 
low-density outer envelope. In contrast, a higher iron 
content in the mantle would increase the density of the 
planet interior, requiring a higher mass fraction of gas. 
While H2 O should dominate any primordial ices forming 
GJ 1214b (contributing more than > 60% by mass), sub- 
stantial quantities (more than 1% by mass) of CO, CO2, 
H 2 S, NH3, N 2 , and CH 3 OH are also expected in ices 
formed fr om a protoplanetary d isk of roughly solar com- 
position (jMarboeuf et alJ l2008). The presence of signifi- 
cant quantities of CO2 could have dramatic effects on the 
evolution and thermal structure of a water p lanet, main- 
tainin g the steam atmosphere in a hot state (|Leger et al.l 
2004). It is important to note that the range of possi- 
bilities for the chemical materials comprising GJ 1214b's 
interior is constrained by the cosmic abundance of the 



elements. We have focused on possible silicate-based in- 
terior composition scenarios for GJ 1214b. Alternatively, 
SiC, graphite, and other carbon compounds could dom- 
inate the interiors of planets formed under condition s 
where C/O > 1 by number (jKuchner k, Seage"rll2005[) . 
but we have not considered this possibility here. 

We find that the choice of opacities used to model 
the GJ 1214b envelope has a considera ble effect on the 
H/He mass fraction constraints in case I. iFreedman et al.l 
(2008) opacities were used to generat e the composi- 
tion c onstraints shown in Figure |3] The Fer guson et al.l 
(2005) opacities, which include condensate and grain 
opacity sources , tend to be higher overall than the 
IFreedman et al.l (|2008l ) opacities, which do not include 
grain opacity sou rces. Consequently, gas envelopes 
modeled with the iFerguson et al.l ([2005D opacities sys- 
tematically require less H/He mass to reproduce the 
observed transit radiu s than those modeled with the 
IFreedman et al.l (|2008|) opacities. This illustrates that 
the relationship between the radial thickness of an H/Hc 
envelope and its mass is sensitive to the precise metal- 
licity and composition assumed. Despite the added un- 
certainty introduced by the metallicity of the H/He layer 
our inferences regarding the plausibility of a primordial 
gas envelope on GJ 1214b remain unchanged. 

6. SUMMARY AND CONCLUSIONS 

The MEarth transiting planet, GJ 1214b, is exciting 
because it lies in a mass and density regime for which 
there are no solar system analogs; GJ 1214b is smaller 
than the ice giants Neptune and Uranus, while larger 
than the terrestrial Earth, Venus, and Mars. We em- 
phasize that, based on its measured planetary mass and 
radius alone, we can constrain GJ 1214b's composition 
but we cannot infer its unique true composition. 

GJ 1214b requires a gas envelope to account for its low 
average density so long as the true planet mass and ra- 
dius lie within 2er of their measured values. With interior 
structure models, we explored three possible scenarios 
for the gas layer and concomitant interior of GJ 1214b. 
An important conclusion from this investigation is that, 
under most of the conditions we considered, GJ 1214b 
would not have liquid water. We summarize more de- 
tailed results for each of the three cases below. 

If GJ 1214b's gas layer was accreted directly from the 
protoplanetary nebula, the primordial H/He layer sur- 
rounding an interior of iron, silicates, and ice would need 
to contain between 0.01% and 5% of the planet mass in 
order to account for the transit radius. This is inter- 
esting because the gas envelope would be less massive 
than Uranus' and Neptune's envelopes (which account 
for 5%-15% of the planet mass), yet greater than Earth's 
or Venus' atmospheres (which contribute 0.0001% and 
0.01% of the planet mass, respectively). 

If, instead, sublimated ices dominate the gas layer, a 
massive water envelope comprising at least 47% of the 
planet mass could account for GJ 1214b's observed pa- 
rameters to within 1 a. We thus do not require an 
H/Hc layer to explain the measured mass and radius. 
In this sublimated ice-dominated case, for our assump- 
tions about the planet albedo (A < 0.35) and internal 
heat flux, we find that GJ 1214b's water envelope would 
generically be too hot to allow liquid water, even at high 
pressures. Instead of a liquid water layer the planet 
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would have a super-fluid water layer sandwiched between 
plasma below and vapor above. To obtain liquid water 
in our model interior, we must decrease GJ 1214b's as- 
sumed intrinsic energy flux to < 0.01 of Earth's intrinsic 
energy flux. 

Third, if the nascent GJ 1214b did not manage to re- 
tain any primordial gas or ices, outgassing from rocky 
material could produce a gas layer surrounding a terres- 
trial interior. In order for sufficient gas to be released in 
this scenario to account for the transit radius, the out- 
gassed atmosphere would need to be hydrogen rich. This 
in turn constrains the mineralogy of the primordial rocky 
material from which GJ 1214b would need to have formed 
in this case. Based on models of outgassing of chondritic 
material found in the solar system, even if hydrogen is the 
dominant outgassed species other heavier volatiles would 
be present as well. It is also expected that atmospheric 
escape would have eroded a substantial amount of the H 
atmosphere (|Charbonneau et al.l 120091) . which provides 
added tension with this scenario. Although outgassing is 
not ruled out as the primary contributor of GJ 1214b's 
gas layer, more modeling is required to verify its viability. 

The ideal hope is that we may gain insights into planet 
migration and formation if future observations succeed in 
constraining the composition of GJ 1214b's atmosphere 
and then by extension the planet's interior. Despite the 
link between atmospheres and interiors, we caution that 
the solar system planet atmospheres, especially terres- 
trial atmospheres that have undergone substantial evo- 
lution, are divorced from their interiors. Progress toward 



constraining GJ 1214b's interior might thus be challeng- 
ing. Additional complications arise from the fact that 
multiple processes may have contributed to the gas layer 
on GJ 1214b. Nevertheless, here are some hopeful pos- 
sible outcomes. If GJ 1214b's gas layer is found to have 
close to solar abundance of He, then it is likely composed 
of nebular gases, and together with atmospheric escape 
estimates may yield bounds on the accretion of gas from 
protoplanetary disks by small planetesimals. If instead 
sublimated ices dominate with no observable H and He, 
it would indicate that GJ 1214b formed beyond its star's 
snow line and migrated inward to its current orbital dis- 
tance. Finally, if GJ 1214b's envelope has no He but 
a substantial amount of H and other volatiles, it would 
likely be the result of outgassing, and may help to con- 
strain its atmospheric mass-loss rate and the broad-brush 
oxidation properties of the rocky material from which 
GJ 1214b formed. No matter which atmospheric sources 
turn out to be dominant on GJ 1214b (be it nebular gas, 
sublimated vapors, or outgassing) we will learn some- 
thing interesting, and perhaps achieve our first glimpse 
of a planetary atmosphere on a whole new class of low- 
mass exoplanets. 
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